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■^J- ' Abstract 

(N ! 

fN| ■ We consider the bosonic sector of a J\f = 2 supersymmetric Chern-Simons-Higgs 

theory in 2 + 1 dimensions. The gauge group is U(l) x SU(N) and has Nf flavors 
of fundamental matter fields. The model supports non-Abelian (axially symmetric) 
vortices when Nf > N, which have internal (orientational) moduli. When Nf > N, 
the solutions acquire additional collective coordinates parameterizing their transverse 
size. We solve the BPS equations numerically and obtain local (Nf = N) and semi-local 
(Nf > N) string solutions. 
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1 Introduction 

Chern-Simons (CS) theories are relevant in a quantum field theory context since they pro- 
vide an alternative gauge-invariant procedure of mass generation pp. Moreover, the high- 
temperature limit of quantum field theories in d = 4 dimensions are effectively three di- 
mensional and CS terms are precisely induced by fermions in d = 3 dimensions through 
the parity anomaly [2j. CS actions play also a role in the analysis of interesting condensed 
matter phenomena [3]- [5], the computation of topological invariants of 3-manifolds [6 J and 
they are connected with d = 2 conformal field theories [7j. 
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CS-Higgs models differ drastically from theories in which solely a Maxwell or Yang- 
Mills term governs the dynamics of the gauge fields. In particular, at the classical level, 
axially symmetric (vortex or string) solutions to the equations of motion necessarily carry 
electric charge [8]-[l5] which, in the non-Abelian case, is quantized (for a complete review 
on CS theories and planar physics see [T6]). 

Bogomolny equations for the non-Abelian CS-Higgs system were first obtained in |15j . 
where it was shown that, as in the Abelian case [12j-[13j, a sixth order potential has to be 
considered. Explicit vortex solutions were exhibited in [15J, with the flux directed in the 
Cartan subalgebra of the non-Abelian group. It is the purpose of the present work to find 
genuine non-Abelian vortex configurations by proposing a more general ansatz as the one 
already considered for Yang-Mills-matter theories p2]-[l9] with gauge group U(l) x SU(N) 
and Nf > N flavors of fundamental matter multiplets. The resulting vortex configurations 
can be characterized by non-Abelian collective coordinates related to orientational degrees 
of freedom and, when Nf > N, to infinitesimal variations of the transverse size. 



2 Model and Notation 

We consider the truncated bosonic sector of the M = 2 SUSY U(l) x SU(N) Chern-Simons- 
Higgs action in 2 + 1 dimensions (the complete SUSY Lagrangian can be found in [21]) 

S = [d 3 x !^ e ^F°A° + ( FLAi - -f^A'AiAf) + 

I ]2 ^ p 2 \ p v p ^ J p v p j 



where e 012 = 1, g 00 = 1, f IJK are the structure constants of the non-abelian group, and the 
quantization condition implies k 2 = m/8ir. Our theory has a sixth-order potential which, 
forced by M = 2 supersymmetry, allows Bogomolny completion of the energy functional with 
all coupling constants at the Bogomolny point (see ref.[22] and references therein for details 
on this point) 
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Here, fj,,v,p = 0, 1, 2 are Lorentz indices, I, J, K — 1, ... , N 2 — l are the SU(N) "color" group 
indices and tj are the anti-hermitian generators of SU(N). The complex scalar multiplets <p{, 
besides the color index i,j,k = 1, . . . , N, possess additional flavor index f,g,h = l,...,Nf 
with Nf > N, thus can be written as iV x Nf matrices. The covariant derivatives and field 
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strengths are defined as 



A SU{N) _ A l 



Tl 



A 4 

Ku = . F i = d \X] + f IJK K A * ■ 

Up to gauge transformations, minima of the potential are given by 

symmetric phase = 

asymmetric phase $ 4>\ = £ diag{l, . . . , 1} . 

In what follows we set, without loss of generality, £ = 1. The energy density is 
and the Euler-Lagrange equations of motion of the theory are 



(3) 

(4) 
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Defining D e = D\ + ieZ)2 with e = ±, and using Gauss' law, we can write the energy as a 
sum of squares 



H 



x 



d 2 x 
D (f) f - ie 
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^ -N]<f>f- ±-4yvT*tf 



+ (D^y (D_ e 4>f)+eV2NF?A 



(9) 



leading to the Bogomolny equations. 



3 Non-abelian local strings 

In this section we investigate the so-called "local Z N string-type solutions" as discussed for 
the Yang-Mills case in [Hj-[l9] (see also [2D]). 

We set Nf = N so the matter fields can be arranged as a square matrix $. The 
Lagrangian (CQ) is then invariant under SU(N) co i or x SU(N) fi avor rotations, 



(10) 
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with U e U(N) loca i and V G SU(N) global . 

We start from the trivial vacuum in the asymmetric phase 

A vac = Q ^ $^ = diag{l,...,l} . (11) 

After a i7(l) x SU(N) gauge transformation we obtain a singular vortex configuration which, 
for the scalar field, takes the form 

^vac y $ = exp ^0 + p T N 2 -l j $ mc = diag { 1; . . . ^ (1 2 ) 

with 

r° = diag{l, . . . , 1} , r^ 2 - 1 = - = diag{l, ...,1,1- N} (13) 

a = -y^2/N erup , 13 = y/2(N-l)/N en V 

Concerning the gauge fields we have 



7Lf c/ ^7^ = 0, A$ Kl ^>Al = 0, I = 1,2,..., TV 2 -2 (14) 

In such configuration, the center of the gauge group SU(N) has been combined with {7(1) 
elements to get a topologically stable string solution possessing both windings, in SU(N) 
and in {7(1) (since tci(SU(N) x U(1)/Z n ) ^ 0, the topology is nontrivial). These kind 
of topological objects are also called Ll Z N strings". In the present case the configurations 
f|T2|) . ffT4"l) represent a (0, . . . , 0, n) singular string. A general (n 1; . . . , n N ) vortex configuration 
can be obtained following the same method. 

This suggests the following ansatz for the regular vortex configuration 

$ = diag{0(r), . . . , 0(r), e"^^(r)} , (15) 



= 7y=/o(r) , 4° = - V + • (16) 

It should be noted that, contrary to the Yang-Mills-Higgs case [19J, the A fields are non- 
trivial. The boundary conditions for the fields are 

<M0) = 0, /(0) = / iv2 - 1 (0) = -en, 

0(oo) = <M°o) = 1, / (oo) = fo N2 -\oo) = f(oo) = f N2 -\oo) = , (17) 
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where we required the solution to be singled- valued at the origin, and to be a pure gauge at 
infinity. With this ansatz both the flux and the energy of these solutions are quantized 
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E = 2nn. 



(18) 



Note that with our conventions n is always a positive integer, and e determines whether the 
flux is positive (e = +), or negative (e = — ). 

Ansatz (1151) corresponds, for the n = 1 case, to what is called an "elementary string" . 
Composite strings can be constructed by introducing windings in several diagonal elements 
in the scalar field and can be seen as the superposition of elementary ones. 

Substituting ansatz (Tl5l)-(fT6l) into the Bogomolny and Gauss equations, we obtain the 
following system of non-linear first-order differential equations 
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(19) 
(20) 
(21) 
(22) 



Note that equations fl22|) can be used to eliminate f and Jq 2 ' 1 in (|20|) and (I2T]) . Equations 
(lT9~j) - (l2"Tj) then correspond to the Bogomolny equations written in the standard form. In 
terms of the profile functions, the energy takes the form 



E = 2tc / r dr 



+ 



(N-l) 

8N 2 k 2 2 



a \2 



((iV-l)0 2 + ^-iV) 2 ((iV 



1)0 2 + fa) 



(tf-faY(^ 2 + (N-l)fa) + 



(N-l) 



AN 2 KiK 2 

((N - 1)0 2 + fa - N) + (N - l)(<9 r 0) 2 + (d r <p N ) 2 
, (N-l) 



(f-f 



N 2 -l 



+ 



(23) 



The magnetic field (which is a pseudoscalar in 2+1 dimensions) and the electric field (with 
component only in the radial direction) take the form 
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E° 



(K) 2 + W) 2 = 7^d r f . 



(24) 



Let us study the solutions of equations (|T9l-(l22l). We can distinguish two cases: first, when 
the U(l) and SU(N) coupling constants are equal and then the effective symmetry group is 
U(N), and second, when the U(l) and SU(N) coupling constants are different. 
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Figure 1: Plot of the magnetic field B° (solid line), electric field E° (dashed line), and energy 
density TL (dotted line) for local vortices with K1/K2 = 1. The distance between the energy density 
maximum and the origin increases with the winding number n. 

3.1 U(N) gauge x SU(N) g i oba i solutions 

If the U{1) and SU(N) coupling constants are equal Ki = k 2 = k, the system of equations 
(H9D-fl22D decouples into 

rd r <f> = -eg<f> 1$.,, = _ * ^ _ i) , (25) 

r 8«r 

rd^N = -eg N2 ~^ N -d^ 1 = - 1) , (26) 

r ok z 

where we have defined 

In terms of the new functions, the boundary conditions are 

0(0) = 0, S(oc)=0, g N *- 1 (0) = -en, ^(oo) = (28) 

0(0) = C, 0(oo) = l, 0tv(O) = O, 0^(oo) = l. (29) 

Since the field has no winding, it does not necessarily vanish at the origin so C is an 
arbitrary parameter. Systems (|25p and (|2"6"|) are formally the same but the functions ((f), g) 
and (4>N,g N _1 ) obey different boundary conditions. Remarkably, both systems coincide with 
those arising in the abelian case discussed in [13] . For the pair (0, <?), boundary conditions 
imply that 

= 1,0 = 0. (30) 

The system ( 126]) for the functions (0jv, g^' 1 ) was solved numerically in [13]. We remark 
that these equations do not depend on N. At the origin, 0jy approaches to zero as r n (this 
fact will be important for the semi-local vortex). 

We show in figure [T] profiles of B°, E° and the energy density 7~C{r). As usually 
happens in CS theories, the sixth order potential makes the maximum of the magnetic field 
to be away from the origin. 
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Figure 2: Profile functions for local vortices, for negative magnetic flux and «i/«2 = 1/2- 

3.2 U(l) gauge x SU(N) gauge x SU(N) g i oba i solutions 

When K\ 7^ k 2 , the complete set of equations (|T9|) - (!22|) has to be solved numerically. We 
used a relaxation method to find explicit numerical solutions. The ratio of coupling constants 
k = k\/k2 is in fact the only independent parameter of the theory, since K\ or k 2 can be 
absorbed by a rescaling. In fact the energy can be expressed in terms of k only. We observe 
that as k departs from 1, / and f N _1 tend to separate from each other as well, forcing 
to be non-constant. As k goes from k < 1 to k > 1, the difference / — f N _1 changes sign, 
forcing the derivative of <p to change sign. Qualitatively, the behavior resembles the U(N) 
case (in which both coupling constant coincide). When varying the winding number n of 
the vortex, the profile functions change in a similar way as they do in the U(N) case. For 
equal n and different N, solutions do not change considerably. We present in figure [2] some 
solutions for the case k = 1/2. 

4 Non-abelian semi-local strings 

We have discussed non-Abelian vortex solutions in a U{l) gauge x SU(N) gauge x SU(N)fi avor 
theory which are usually called local vortex solutions. In order to have semilocal vortex 
solutions, those for which Nf > N, one has to extend the matter content of the theory by 
adding N e extra flavors so that Nf = N + N e [23]- [21] . For definiteness we take N e equal 
to iV (but a general case can be equally treated). Then, the symmetry of the model is 
U(l) gauge x SU(N) gauge x SU(2N)fi avor . We call X the extra matter fields. 
In this case the trivial vacuum state is given by 

$ mc = diag{l, . . . , 1} , X vac = diag{0, . . . , 0} , A™ c = , (31) 

and it is invariant under the following transformation 

QVOC _^ y-lQvacy ^ ^vac _^ y-l^vacy = q ^ jyac _^ y-l^acy = g ^ ^ 

where we have chosen, as in the local case, a gauge element U = V^ 1 with V a iV x N flavor 
block and we have called V the N x N flavor block acting on the extra X mater fields. 
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We follow now the same steps as in the local case: first we perform a rotation of the 
vacuum to find a singular vortex configuration and then propose an ansatz for a regular 
configuration which we write explicitly 

$ = diag{0(r), . . . , 0(r), e"^ ne <Mr)} , X = diag{x(r), . . . , x(r), X7v(r)} , 

A * = jf £ ^ ne + /^M) diag{l, 1, . . . , 1 - N} - ± £ ^(ne + /(r)) Id , 

Ao = ^ ff~\r) diag{l, -N} + ^ / (r) Id (33) 

where Id is the N x N identity matrix. It should be stressed that more general ansatz 
could lead to solutions which exhaust the number N p of collective coordinates related to the 
transverse moduli space [21]. Our ansatz corresponds to just two collective coordinates. 

Substituting (1331) into the Bogomolny equations, we arrive to a system of six differ- 
ential equations for the fields 0, 0at, x, xn, f and f N ~ 1 , and two constraints for the fields 
/o and /q -1 . The solutions to these equations will have the same energy and flux than the 
local ones (ITS]) , provided they satisfy the same boundary conditions (TT71) . At infinity the 
solutions should reach the vacuum state, so we impose 

0(oo) = 0jv(oo) = 1, x (oo) = xw(oo) = , 

/o(oo) = /o^-^oo) = /(oo) = /^(oo) = . (34) 

Again we distinguish the case in which the U(l) and SU(N) gauge coupling constants 
are equal from the one in which they are different. 



4.1 U(N) gauge x SU(2N) flavor semilocal solutions 

As it happens for local vortex, when Ki = k 2 = k the equations of motion decouple into two 
sets of independent equations. One corresponds to a system with no winding 

rd r (f) = -eg<fi , rd r x = -egx 

ld r g = -^- 2 (<p 2 + x 2 )(<P 2 + X 2 -l) , (35) 
and the other corresponds to a system with winding n 

rd r (f) N = -eg N2 ~ l (f) N , rd r XN = -e (g 1 * 2 ' 1 + en J xn 

-d r g N2 - 1 = ~ + X%) (<t> 2 N + X%-1) (36) 



(we have defined g and g N _1 as in (J2] 

As in the case of local strings, these equations coincide with those that arise in the 
abelian case |25j. The system (135]) admits a trivial solution 

g = x = , = 1, (37) 
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Figure 3: The magnetic field for semi- local vortices with k% = k,2 for different p values 



while (I36I) can be solved numerically. Combining the equations for \n and 4>n, we get 

Xn = P — , (3t 
where p is a (complex) integration constant. So, we finally have 



rd r (j) N = -eg 



N 2 -l 



V"- 1 = Cfl + A - l") (l + 4 ) A • (39) 



The same argument applied to \ and <fi determines that \ = a<f), but since <f)(oo) = 1 and 
x(oo) = 0, a must vanish, and so % = everywhere. The solutions we obtain are then the 
most general ones. The energy is independent of the complex parameter p, which is then 
associated with the "size moduli". One has has in general, Hi complex parameters p« for 
arbitrary flavor Nf > N . For elementary strings one can see that n« = Nf — N (Since we 
are considering Nf = N + 1, we have in this case just one complex parameter). 

From figure [3] we see that the solutions spread when the parameter \p\ is incremented. 
As p increases, since the flux is conserved, the extremum of the magnetic field must approach 
to the origin to compensate the spread towards spatial infinity. A similar phenomenon occurs 
with the energy. 



4.2 U (l)gau 9 e x SU (N) gauge x SU(N + 1) fi avor semilocal solutions 

We consider now the case in which K\ ^ /? 2 . The same arguments as above lead to x{ r ) — 0, 
so that the flavor group is in fact reduced from SU(2N) to SU (N + l). Similarly, the relation 
(138]) holds, and the differential equations do not decouple. The system is very analogous to 
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the system (|T9l) - (l22l) obtained for the local case, with the only difference that, in the equations 
for the /'s, the field 4>n gets locally scaled 



Then, semi-local solutions arising from these equations are similar to those shown in figure [2] 
for the local case, with the only difference that they are smoother since they decay as powers 
of r at spatial infinity. 

In summary, our main task in this work was to solve the BPS equations for a non- 
Abelian Chern-Simons-Higgs theory. By proposing an axially symmetric ansatz we obtained 
non-Abelian vortex solutions and discussed their properties. A class of vortex solutions in 
non-Abelian CS theories were already known [9J-[10],[TT], [IS]. The solutions discussed here 
correspond to more fundamental vortices in the sense that they are genuinely non-Abelian 
while the former correspond to vortices with the gauge flux in the Cartan algebra of 



The model discussed here is indeed related to the one analyzed by [H]-[l9] except 
that in our case the dynamics of the gauge fields is governed by a CS action instead of a 
Yang-Mills one. This drastically changes the vortex properties, in particular forcing them to 
carry electric charge. 

In the case of local vortices, our solutions generalize those discussed in |20j to the case 
in which the gauge group is U(l) x SU (N), with the Abelian and non-Abelian sectors having 
different gauge coupling constants. When both couplings are equal, the equations decouple 
into two sets of equations that coincide with those arising in the Abelian case. This is not the 
case when the couplings are different and the BPS equations do not decouple. Nevertheless 
we were able to construct explicit solutions and discuss their properties. 

Further, we have also considered semi-local vortices, by allowing the flavor number Nf 
to be larger than the color number iV c . As already noted for the Yang-Mills-Higgs system, the 
main feature in this case is that the solutions develop an additional moduli p related to the 
vortex transverse size, thus modifying the asymptotic behavior of the fields, from exponential 
decay for the case of local vortices, to a decay as negative power of the radial coordinate for 
the semi-local ones. Interestingly enough, one can see from our explicit solutions how the 
size of the vortex grows with p. 
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